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The optimal reactive power dispatch is a kind of optimization problem that 
plays a very important role in the operation and control of the power system. 
This work presents a meta-heuristic based approach to solve the optimal 
reactive power dispatch problem. The proposed approach employs Crow 
Search algorithm to find the values for optimal setting of optimal reactive 
power dispatch control variables. The proposed way of approach is 
scrutinized and further being tested on the standard IEEE 30-bus, 57-bus and 
118-bus test system with different objectives which includes the 
minimization of real power losses, total voltage deviation and also the 
enhancement of voltage stability. The simulation results procured thus 
indicates the supremacy of the proposed approach over the other approaches 


dispatch cited in the literature. 
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1. INTRODUCTION 

Optimal Reactive Power Dispatch (ORPD) enhances system security and helps in proper planning 
of reactive power and its dispatch in the system. The ORPD problem, which is a non-linear multi-objective 
optimization problem aims at minimizing the objectives such as real power loss, total voltage deviation and 
enhancing voltage stability by optimal setting of the control variables such as voltages of generator, 
transformer Tap set values and the reactive power output obtained from shunt capacitors in an optimal way 
which has to meet the required equality and inequality constraints. Over the last few decades, the ORPD 
problems were solved using various traditional methods like Newton's approach, linear programming, 
interior point method, etc. [1]-[4]. There are certain disadvantages also pertaining to the traditional 
methodologies, which are difficulties in finding a global optimal solution, insecure convergence, complexity 
and bad initial termination. 

In order to overcome the above illustrated drawbacks, heuristic methodologies have been under 
research for solving ORPD problem. In the past few decades, heuristic algorithms such as Genetic Algorithm 
[5]-[9], Evolutionary Programming [10], [11], Particle Swarm Optimization [12]-[15], Ant Colony 
Optimization [16], [17], Bacterial Foraging Optimization [18], Seeker Optimization Algorithm [19], 
Differential Evolution [20], Gravitational Search Algorithm [21], Cuckoo Search Algorithm [22] were used 
for solving ORPD problem. The above stated heuristic algorithms have overcome the drawbacks in the 
traditional methods, but also have certain limitations, that is they easily get trapped in the local optima and 
premature convergence would occur . In this paper, a CSA [23] based approach is being proposed to solve 
the ORPD problem. The problem is formulated as a nonlinear optimization problem with both equality and 
inequality constraints. In this study, different objectives are considered such as minimizing the power losses, 
improving the voltage profile and enhancing power system voltage stability. The proposed way of approach 


Journal homepage: http://iaescore.com/journals/index.php/IJECE 


1424 O ISSN: 2088-8708 


has been examined and tested on the standard IEEE 30-bus, 57-bus and 118-bus test system [24]-[26]. The 
prospective and effectiveness of the proposed approach is demonstrated and the results are compared with 
those cited in the literature. 


2. ORPD PROBLEM FORMULATION 
In the present work the ORPD problem is formulated as a multi-objective optimization problem 
which is listed below: 


2.1. Minimization of real power loss 
The objective of ORPD problem is minimization of real power loss while satisfying its various 
equality and inequality constraints. 


N, 
Minimize F, =) g8,[(V? +V? -2V,V, cos(6, -6,))] 
al (1) 


where N; represents the number of transmission lines, g,is the conductance of i™ transmission line, V; and V; 
is the voltage magnitudes of i" and j" buses, 5 and 5 is the voltage phase angle of i" and j” buses. The state 
r F 


variables for ORPD problem are, 


AATA 


ETERA 2) 
where Pg; denotes the slack bus power, V; denotes the load bus voltages, Qg denote the reactive power 
output of the generators, Ng denotes the number of voltage controlled buses, Npg denotes the number of load 
buses. The control variables for ORPD problem are, 


i VGN One O get Te 3) 


where Nr and Nc denotes the number of tap changing transformers and shunt VAR compensators. Vg 
denotes the voltages of voltage controlled bus and T denotes the transformer tap ratio and Qc denotes the 
reactive power output of shunt VAR compensators. 


2.2. Minimization of total voltage deviation 
The minimization of voltage deviation improves the voltage profile of the system, thereby 
enhancing the security and quality of the system. The objective function is given as, 


Minimize F,= > |V,-V/“ | (4) 
i€Npo 
where V; is the reference value of i™ bus voltage magnitude which is usually 1.0p.u. 


2.3. Enhancement of voltage stability 

Voltage stability is the ability of the system to maintain its voltages within its permissible limits. 
Voltage instability occurs only when a system is subjected to disturbances in the system. Voltage stability 
can be improvised by minimizing the voltage stability indicator L-index at all buses. L-index is usually in the 
range of 0 to 1 for all load buses. The L-index at a bus denotes the chances of the voltage collapse condition 
of that bus. L-index L; of the j” bus is given as, 


7 ST x Where j = 1,..., Neo (5) 
i=l j 
Fi 7 ly, T [r] (6) 


where Npy denotes number of PV buses. Y; and Y, are the sub-matrices that are obtained after separating PQ 
and PV bus parameters, 
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L-indéx is calculated for all PQ buses. The L-index value can be described as, 


(7) 


Lx = Max(L,), where j = J,..., Npg (8) 


max 


A lower value of L denotes a stable system. In order to improve the voltage stability and to move the 
system from voltage collapse point, the objective can be described as, 


Minimize F, = Lax 


(9) 
where Lmax is the maximum value of L-index. 
2.4. Constraints 
The objective functions are subject to both Equality and Inequality constraints as below: 
2.4.1. Equality constraints 
Nb 
P; — Po, — XV V,[G, cos(6, -5,) + B, sin(6, -5,)]=0 i=1,2,---, Nb 
jh (10) 
Nb 
06; -— Qn; — X VV [G; sin(6, -8,)- B, cos(6, -6,)]=0 i=1,2,---,Nb 
j=l (11) 


where Nz is the number of buses, Pg; and Qa; are generated active and reactive power, Pp; and Qp; are active 
and reactive power of load, Gj is the transfer conductance and B; is the transfer susceptance between i and 


2.4.2. Generator constraints 
Generator voltages and reactive power output have to be within its permissible limits described as 
follows, 


Van" SVa svan i= 1. Ng (12) 


On" < Qo SORS, i=1,...,Ng (13) 


2.4.3. Transformer constraints 
The transformer tap settings have to be within its lower and upper limits as follows, 


e EGS Dis i=1,..., Nr (14) 


2.4.4. Shunt VAR compensator constraints 
Shunt VAR compensators should be within its lower and upper limits as follows, 


a Oa On i=1,.., Ne (15) 


2.4.5. Security constraints 
The load bus voltages and transmission line loadings have to be within its prescribed limits, 


yR" <V; av i= 1,...,Npo (16) 
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Thus the generalized objective function can be formulated as, 


NPQ - NG = BE Nir holed 
Min Obj= Fajt A, Xv, Vey + dy Y (Qe: 7 Qu)’ Ak; Ae (Si: -Su 
i=l i=l i=l (18) 


where Fos; is the objective function, Ay, A, and Ag are the penalty factors. 


3. CROW SEARCH ALGORITHM 

Crows or Corvids are intellectual omnivores; natural history books remain to be an evidence for it. 
Crows have remarkable abilities like problem-solving skills, communication skills and adaptability. Crows 
are known for its excellence in memory, certain vital researches show that crows don’t forget a face and 
hence alerts other crows how to identify the individuals. Certain behavior of crows is enlisted [23], 
a. Crows live in groups 
b. Crows have excellent memory on their position of hidden places 
c. Crows follows each other to perform acts of thievishness 
d. Crows hide their collectives that have been theft 

Crow Search Algorithm (CSA) is developed based on the above nature and behavior of crows. The 
algorithm has d-dimensional environment with N number of crows and the position of crows ( y+) which can 


be specified by a vector, 
MOS eee | (19) 


where i = 1,2,...,N ; k =1,2....,iter™* ; iter™ is the maximum number of iterations 


In accordance with its memory capacity, the algorithm proceeds as, at k" iteration, the position of 
hiding place of i™ crow is given by, M;. For better illustration, assume that j crow wants to visit its hiding 
place at k” iteration, at this instant of iteration, i" crow follows j™ crow to know its hidden place, here there 
are two possibilities, 

Possibility 1: The crow j being unaware of crow i, shows its hidden place, hence at this instant the new 
position of crow i become, 


xe = $i 4 r, x fli’ x(M i" — xi) (20) 


where r; is a random number with uniform distribution between 0 and 1, fl;"*" denotes the flight length of 
crow i at iteration iter. The value of fl has great impact on the search space of the algorithm, if f! is a smaller 
value than it results in local search and if fl is a larger value it results in global search. 

Possibility 2: The crow j aware of crow i that it is following it, in order, hence to protect its collect from crow 
i, crow j will move to another position to divert crow j, then the new position is thus given by, 


te iter iter it iter 

X eril = Xe +r x fle x(M 5" =A )r, 2 AP” (21) 
! a random position otherwise 

where AP}! denotes the awareness probability of crow j at iteration iter. This factor decides whether the 

search space is intensified or diversified. When AP is increased, the search space gets increased thereby, 

results in global optimal and vice versa. 


4. APPLICATION OF CROW SEARCH ALGORITHM FOR ORPD PROBLEM 
The sequence of steps that ought to be followed in the implementation of the CSA is given in this 
section. 
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Step 1: Initialization of algorithm parameters and constraints 

The algorithm parameters comprises of population size (N), maximum number of iterations (iter), 
flight length (fI) and awareness probability (AP) and the constraints include power balance equality 
constraints, line flow and voltage constraints. 
Step 2: Initialization of the position and memory of crows 

The N population of crows is randomly positioned in a d-dimensional search space. Each crow 
denotes a possibility of feasible solution of the problem and d is the number of control variables which 
includes generator voltages, transformer tap settings and reactive power output of shunt capacitor. The 
memory of each crow is initialized. At the beginning of iteration iter, it is assumed that the crows have 
hidden their foods at their initial positions. 
Step 3: Evaluate fitness (objective) function 

For each crow, the position is determined by fitting the control variable values into the objective 
function (minimization of real power loss, total voltage deviation and voltage stability indicator). 
Step 4: Generate new position 

Crows finds a new position in the d-dimensional search space by as follows: suppose crow i wants 
to find a new position. For this, the crow randomly selects one of the crows, let that be crow j and follows it 
to discover the Position of collecting hidden by this crow (m). The new position of crow i is given by 
Equations (20) and (21). 
Step 5: Check the feasibility of new positions 

The viability of the new position of each crow thus obtained is checked and the position is updated 
based on it. If the new position, thus obtained is not viable, then the crowd stays in the current position and 
does not move to the new position found. 
Step 6: Evaluate the fitness function of new positions 

The fitness function i.e. objective function value for the new position of each crow is evaluated. 
Step 7: Update memory 

The crows update their memory as follows: 


sisi oe is better than fm’) (22) 
m’ = ae 
liter 


m otherwise 


where fo; denotes the objective function value. 

It is seen that if the fitness function value of the new position of a crow is better than the fitness 
function value of the memorized position, the crow updates its memory by the new position. 
Step 8: Check termination criterion 

Steps 4 to 7 are repeated until maximum iteration is reached. When the termination criterion is met, 
the best position of the memory in terms of the objective function value is reported as the solution of the 
optimization problem. 


5. RESULTS AND DISCUSSION 

The present work is being tested on standard IEEE-30, 57 and 118 bus systems and the results are 
obtained. The description of these studied test systems is depicted below. The software is written in 
MATLAB R2015 computing environment. The various algorithm parameters are initialized and are set to be 
as: The value of flock size (population) is set to 75, the awareness probability index determines whether the 
search space is intensified or diversified and is set to 0.5, the flight length is assumed to be 2 and the 
maximum number of iterations performed is set to 200 for all the test cases considered. The results of interest 
are boldfaced in the respective tables to indicate the optimization capability of the proposed algorithm. 


5.1. Case-1: Minimization real power loss 

In this case, the proposed algorithm is executed considering the minimization of real power loss 
alone as the objective function. The convergence characteristic of the algorithm considering the real power 
loss is shown in Figure 1, which indicates fast and smooth convergence of CSA. The superiority of the 
aforesaid CSA based approach for solving ORPD problem can be witnessed from the comparison made 
between other optimization techniques from Table 1, Table 2 and Table 3. The best power loss obtained 
using CSA for IEEE 30, 57 and 118 bus systems are 2.8507 MW, 15.1934 MW and 76.7783 MW 
respectively, which is lesser than result reported in [12], [14], [19], [21]. 
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Table 1. Comparison of results for minimization of active power loss for IEEE-30 bus system 
Methodology CSA CLPSO [14] GSA [21] 
Power loss (MW) 2.8507 4.5615 4.5143 


Table 2. Comparison of results for IEEE 57-bus system 
Methodology CSA CLPSO [14] SOA [19] GSA [21] 
Power loss (MW) 15.1934 24.5152 24.2654 23.4611 


Table 3. Comparison of results for IEEE 118-bus system 
Methodology CSA PSO [12] SOA [19] GSA [21] 
Power loss (MW) 76.7783 131.99 114.9501 127.7603 


5.57 | 4 


45h | 4 


Power Loss (MW) 
a 
T 


25 r r r r r r r r r 
0 20 40 60 80 100 120 140 160 180 200 
Iteration 


Figure 1. Convergence characteristics considering the real power loss as objective 


5.2. Case-2: Minimization of total voltage deviation 

The proposed CSA approach is also applied for minimization of total voltage deviation of IEEE-30 
bus test network and the result yielded from this approach is illustrated in Table 4 and are compared with 
those reported in the literature. The minimum total voltage deviation obtained by the proposed method is 
0.0907, which is lesser than results reported in [12], [14]. The convergence characteristic of voltage deviation 
versus number of iterations is depicted in Figure 2. 


0.4 


Voltage Deviation 


r r t r r 
20 40 60 80 100 120 140 160 180 200 
Iteration 


Figure 2. Convergence characteristics considering voltage deviation as objective 


Table 4. Comparisons of results for voltage profile improvement IEEE-30 bus system 
Methodology CSA PSO [12] CLPSO [14] 
x Voltage deviation (p.u) 0.0907 0.2450 0.2577 
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5.3. Case-3: Enhancement of voltage stability 

In this case the enhancement of voltage Stability is taken as objective function. The solution 
obtained by the proposed method and reported in literature by other methods is illustrated in Table 5. The 
voltage stability indicator obtained by the CSA is 0.1180, which is lesser than results reported in [20], [21] 
and which is proving the excellence of the aforesaid CSA algorithm over other optimization techniques. The 
convergence characteristic of L-index versus number of iterations are depicted in Figure 3, which shows fast 
and smooth convergence characteristics of CSA. 


Table 5. Comparison of results for enhancement of voltage stability IEEE-30 bus system 
Methodology CSA DE [20] GSA [21] 
L-index (p.u) 0.1180 0.1246 0.1368 


0.134 ' r r r r r - - - 
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Figure 3. Convergence characteristics considering voltage stability indicator as objective 


The optimal setting of the control variables for the IEEE-30 bus system is illustrated in Table 6. The 
optimal control variable setting for IEEE-57 and IEEE-118 bus system for case-1 (minimization of real 
power loss) is illustrated in Table 7 and Table 8 respectively. 


Table 6. Control variable settings for IEEE-30 Bus System 


Control variable Case-1 Case-2 Case-3 
Vai (p.u) 1.1000 1.0152 1.1000 
Va (p.u) 1.0975 1.0006 1.0882 
Vas (p.u) 1.0796 1.0173 1.1000 
Vas (p.u) 1.0867 1.0027 1.0885 
Vai (p.u) 1.1000 1.0736 1.1000 
Vai3 (p.u) 1.1000 1.0172 1.1000 

Too 1.0665 1.0961 1.0025 
To-10 0.9000 0.9000 0.9000 
Tair 0.9880 0.9972 0.9675 
T2s-27 0.9738 0.9692 0.9078 
Qcio (MVAR) 5.0000 4.0381 5.0000 
Qei2 (MVAR) 5.0000 4.7556 5.0000 
Qcis (MVAR) 5.0000 4.9998 4.3599 
Qci7 (MVAR) 5.0000 0.0006 4.9892 
Qce2 (MVAR) 4.0451 4.9979 4.8982 
Qooi (MVAR) 5.0000 4.9785 0 
Qcz (MVAR) 2.6117 5.0000 0 
Qca (MVAR) 5.0000 5.0000 0 
Qc (MVAR) 2.2796 2.8054 0 

Power loss (MW) 2.8507 10.3406 9.0087 
= Voltage 2.0447 0.0907 2.3063 
deviation 

L-index 0.1261 0.1489 0.1180 
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Table 7. Optimal settings of control variables for IEEE 57-bus system 


Control Control Control Control Control 


. Value . Value : Value f Value : Value 
variables variables variables variables variables 
Voi (p.u) 1.0468 Vo. (p.u) 1.0290 Qc3 (MVAR) 13.8328 Tia 0.9001 Ti1-43 0.9564 
Voz (p.u) 1.0457 Tais 0.9964 To4-25 1.1000 T15-45 0.9634 T40-56 0.9874 
Vo3 (p.u) 1.0423 Tais 0.9708 To4-25 0.9922 T1446 0.9573 T3957 0.9801 
Voce (p.u) 1.0550 T21-20 1.0018 T24-26 1.0281 Tio-51 0.9667 To-55 0.9987 
Vas (p.u) 1.0638 Qci (MVAR) 9.0156 T7-29 0.9924 Ti3-49 0.9239 
Voo(p.u) 1.0346  Qo(MVAR) 17.3349 Tas 0.9705 Power loss (MW) 15.1934 


Table 8. Optimal settings of control variables for IEEE 118-bus system 


Control Value Control Value Control Value Control Value Control Value 
variables variable variables variables variables 


Voi (p.u) 1.0238 Vo34 (p.u) 1.0175 Vo70 (p.u) 1.0259 Vo103 (p.u) 1.0307 Qc (MVAR) 39.3480 
Vas (p.u) 1.0287 Voao (p.u) 1.0098 Von (p.u) 1.0250 Va105 (p.u) 1.0154 Qces3 (MVAR) 8.6446 
Vos (p.u) 1.0764 Voa2 (p.u) 1.0149 Vor (p.u) 1.0146 Vo107 (p.u) 1.0114 Qeios (MVAR) 0.6263 
Vaio (p.u) 1.0946 Voae (p.u) 1.0262 Vor (p.u) 1.0121 Voi 10 (p.u) 1.0181 Qeio7 (MVAR) 26.8462 
Vaiz (p.u) 1.0247 Voao (p.u) 1.0410 Von (p.u) 1.0314 Voi 11 (p.u) 1.0246 Qeiio (MVAR) 13.4358 


Vais (p.u) 1.0207 Vos4 (p.u) 1.0237 Vaso (p.u) 1 0452 Voi 12 (p.u) 1.0145 Ts. 1.0150 
Vais (p.u) 1.0239 Voss (p.u) 1.0249 Vass (p.u) 1.0564 Voi13(p.u) 1.0322 T25-26 1.0500 
Vaig (p.u) 1.0181 Vase (p.u) 1.0234 Vos7 (p.u) 1.0581 Vout6 (p.u) 1.0299 T1730 1.0354 
Voz (p.u) 1 .0474 Vaso (p.u) 1 04 13 Voso (p.u) 1.0722 Qc34 (MVAR) 6.8484 T37-38 1 0 1 37 
Vo2s5 (p.u) 1.0765 Voe1 (p.u) 1.0345 Voo0 (p.u) 1.0479 Qca (MVAR) 2.2374 T59-63 0.9791 
Vazs (p.u) 1.1000 Voo2 (p.u) 1.0327 Voo1 (p.u) 1.0465 Qcas (MVAR) 23.4731 To1-64 0.9991 
Voz (p.u) 1.0345 Vos (p.u) 1.0354 Voo2 (p.u) 1.0550 Qc (MVAR) 0 T65-66 0.9668 
Vo31 (p.u) 1.0221 Veo (p.u) 1.0517 V a99 (p.u) 1 .0368 Qcag (MVAR) 9.6476 T68-69 0.9380 
Vo32 (p.u) 1.0296 Veo (p.u) 1.0546 Vo100 (p.u) 1.0411 Qca (MVAR) 1 1.7550 Tso-81 0.9742 
PL (MW) 76.7783 


6. CONCLUSION 

In this work, the CSA algorithm was proposed for solving optimal reactive power dispatch problem 
and successfully implemented in IEEE 30, 57 and 118 bus systems. The results obtained from the CSA 
approach were compared with those reported in recent literature and hence the CSA algorithm proves its 
capability of solving ORPD more efficiently in terms of its search capability and robustness. The supremacy 
of CSA over the other approaches was observed. In accordance with the results obtained, the CSA algorithm 
has a simple structure and quick convergence characteristics and therefore can be used to solve ORPD in 
large scale power systems and may be recommended as a very promising algorithm for solving complex 
engineering optimization problems. 
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